Integrability Test for Discrete Equations via Generalized Symmetries

In this article we present some integrability conditions for partial difference equations obtained using the formal symmetries approach. We apply them to find integrable partial difference equations contained in a class of equations obtained by the multiple scale analysis of the general multilinear dispersive difference equation defined on the square.Comment: Proceedings of the Symposium in Memoriam Marcos Moshinsk

Probing the dynamics of identified neurons with a data-driven modeling approach

In controlling animal behavior the nervous system has to perform within the operational limits set by the requirements of each specific behavior. The implications for the corresponding range of suitable network, single neuron, and ion channel properties have remained elusive. In this article we approach the question of how well-constrained properties of neuronal systems may be on the neuronal level. We used large data sets of the activity of isolated invertebrate identified cells and built an accurate conductance-based model for this cell type using customized automated parameter estimation techniques. By direct inspection of the data we found that the variability of the neurons is larger when they are isolated from the circuit than when in the intact system. Furthermore, the responses of the neurons to perturbations appear to be more consistent than their autonomous behavior under stationary conditions. In the developed model, the constraints on different parameters that enforce appropriate model dynamics vary widely from some very tightly controlled parameters to others that are almost arbitrary. The model also allows predictions for the effect of blocking selected ionic currents and to prove that the origin of irregular dynamics in the neuron model is proper chaoticity and that this chaoticity is typical in an appropriate sense. Our results indicate that data driven models are useful tools for the in-depth analysis of neuronal dynamics. The better consistency of responses to perturbations, in the real neurons as well as in the model, suggests a paradigm shift away from measuring autonomous dynamics alone towards protocols of controlled perturbations. Our predictions for the impact of channel blockers on the neuronal dynamics and the proof of chaoticity underscore the wide scope of our approach

Lie discrete symmetries of lattice equations

We extend two of the methods previously introduced to find discrete symmetries of differential equations to the case of difference and differential-difference equations. As an example of the application of the methods, we construct the discrete symmetries of the discrete Painlev\'e I equation and of the Toda lattice equation

Classification of five-point differential-difference equations

Using the generalized symmetry method, we carry out, up to autonomous point transformations, the classification of integrable equations of a subclass of the autonomous five-point differential-difference equations. This subclass includes such well-known examples as the Itoh-Narita-Bogoyavlensky and the discrete Sawada-Kotera equations. The resulting list contains 17 equations some of which seem to be new. We have found non-point transformations relating most of the resulting equations among themselves and their generalized symmetries.Comment: 29 page

Marginal cost-based pricing of distribution: a case study

This paper presents results of a software development project carried out by the “Electricity North West” (ENW) and “TNEI” to find economic use-of-system charges for the extra high-voltage (EHV) network. Several cost-based charging models which satisfy principles set by the Regulator, such as cost reflectivity, predictability, stability and transparency were developed. In this paper, the emphasis is put on the developed software and the comparison of nodal marginal charges obtained from the proposed pricing models

QCDF90: Lattice QCD with Fortran 90

We have used Fortran 90 to implement lattice QCD. We have designed a set of machine independent modules that define fields (gauge, fermions, scalars, etc...) and overloaded operators for all possible operations between fields, matrices and numbers. With these modules it is very simple to write high-level efficient programs for QCD simulations. To increase performances our modules also implements assignments that do not require temporaries, and a machine independent precision definition. We have also created a useful compression procedure for storing the lattice configurations, and a parallel implementation of the random generators. We have widely tested our program and modules on several parallel and single processor supercomputers obtaining excellent performances.Comment: LaTeX file, 8 pages, no figures. More information available at: http://hep.bu.edu/~leviar/qcdf90.htm

The Taming of QCD by Fortran 90

We implement lattice QCD using the Fortran 90 language. We have designed machine independent modules that define fields (gauge, fermions, scalars, etc...) and have defined overloaded operators for all possible operations between fields, matrices and numbers. With these modules it is very simple to write QCD programs. We have also created a useful compression standard for storing the lattice configurations, a parallel implementation of the random generators, an assignment that does not require temporaries, and a machine independent precision definition. We have tested our program on parallel and single processor supercomputers obtaining excellent performances.Comment: Talk presented at LATTICE96 (algorithms) 3 pages, no figures, LATEX file with ESPCRC2 style. More information available at: http://hep.bu.edu/~leviar/qcdf90.htm

Coseismic horizontal slip revealed by sheared clastic dikes in the Dead Sea Basin

Peer reviewedPostprin

Tunneling between single and multi-centered black hole configurations

We find a gravitational instanton that connects an initial state corresponding to a single-centered extremal Reissner-Nordstrom (ERN) black hole configuration, to a final state corresponding to a multi-centered configuration. This instanton is interpreted as describing quantum tunneling between the two different black hole solutions. We evaluate the Euclidean action for this instanton and find that the amplitude for the tunneling process is equal to half the difference in entropy between the initial and final configurations.Comment: 8 pages, 4 figures. v4: final version accepted for publication in Phys. Rev.